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Abstract. The recently introduced multifractal network generator (MFNG), has 
been shown to provide a simple and flexible tool for creating random graphs with 
very diverse features. The MFNG is based on multifractal measures embedded in 
2d, leading also to isolated nodes, whose number is relatively low for realistic cases, 
but may become dominant in the limiting case of infinitely large network sizes. Here 
we discuss the relation between this effect and the information dimension for the Id 
projection of the link probability measure (LPM), and argue that the node isolation 
can be avoided by a simple transformation of the LPM based on rotation. 



1. Introduction 

The network approach for describing complex natural, social and technological 
phenomena has become very popular in the recent years. This can be accounted for the 
generality of its fundamental concept of representing the connections among the units 
(building blocks) of the system under study with a graph [H |2] . Over the last decade 
it has turned out that networks corresponding to realistic systems can be highly non- 
trivial, characterised by a low average distance combined with a high average clustering 
coefficient [3], anomalous degree distributions [H [5] and an intricate modular structure 

From the beginning of this new interdisciplinary field, network models have been 
playing a crucial role since they enable singling out the simplest aspects of complex 
structures and, thus, are extremely useful in understanding the underlying principles. 
Furthermore, models can also help testing hypotheses about measured data. In parallel 
with the discovery of the fine structure of real networks, many important and successful 
models have been introduced over the past 10 years for interpreting the different aspects 
of the studied systems. However, most of these models explain only a particular aspect of 
the network (clustering, a given degree distribution, etc.), and for each newly discovered 
feature a new model had to be constructed. 

Due to this proliferation of network models, the concept of general network models 
and methods for generating graphs with desired properties has attracted great interest 
lately. A number of noteworthy methods have been proposed starting from the 
exponential random graph model [9], [101 HH [I2], through the hidden variable models 
[T3| [HI (including the systematic study of the entropy of network ensembles [15] ) , the 
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dK series approach [15] and the use of p-adic randomised Parisi matrices (TTJ HB] to the 
Kronecker-graph approach [191 120] • 

A very recently introduced approach along this line is the multifractal network 
generator (MFNG) [21], which was shown to be capable of generating a wide variety 
of network types with prescribed statistical properties, (e.g., with degree- or clustering 
coefficient distributions of various, very different forms). At the heart of this method 
lies a mapping between 2d measures defined on the unit square and random graphs. 
The main idea is to iterate a suitably chosen self similar multifractal (becoming singular 
in the limiting case) and enlarge the size of the generated graph (becoming infinite in 
the limiting case) in parallel. A very unique feature of this construction is that with the 
increasing system size the generated graphs become topologically more structured. 

However, a slight drawback of the method is that when the size of the generated 
networks (and in parallel, the number of iterations in the multifractal) grow to infinity, 
isolated nodes overtake the majority of the graph in most settings, (see the SI of Ref.[2T] 
for more details). Although this effect is usually not an issue when constructing graphs 
of sizes comparable to real networks, finding a way to circumvent it would still provide a 
noteworthy improvement, especially in the light of the non-trivial connections between 
convergent graph sequences in the infinite network size limit and 2d functions on the 
unit square [221 123] . 

In this article we study the relation between the node isolation effect and the 
Id projection of the link probability measure (LPM) used in the graph generation 
process. Furthermore, we propose a natural method to overcome the problem by 
a simple transformation based on rotation. The paper is organised as follows. In 
Sect 121 we overview the definition and most important properties of the multifractal 
network generator, while in Sect 131 we discuss the connection between node isolation 
and multifractality. We continue by proposing a modification of the original method 
avoiding the node isolation in SectlH, which is tested in practise in SectEl Finally, we 
conclude in SectlHl 

2. The multifractal network generator 

The multifractal network generator was inspired by earlier results from L. Lovasz and co- 
workers proving that in the infinite network size limit, a dense graph's adjacency matrix 
can be well represented by a continuous function W{x,y) on the unit square [221 123] . 
A similar approach was introduced by Bollobas et al. in Refs. [211 125] and was used to 
obtain convergence and phase transition results for inhomogeneous random (including 
sparse) graphs. This two variable symmetric function (which can have a very simple 
form for a variety of interesting graphs, and was supposed to be either continuous or 
almost everywhere continuous) predicts the probability whether two nodes are connected 
or not. 

In case of the MFNG the mentioned W{x, y) is replaced by a self-similar multifractal 
[2T] . We start by defining a generating measure on the unit square by dividing identically 
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both the X and y axis to m (not necessarily equal) intervals, splitting it to rectangles, 
and assigning a probability pij to each rectangle G [1, m] denote the row and column 
indices). The probabilities are assumed to be normalised, — 1 symmetric 

Pij ~ Pji- (Note that we normalize the probabilities of the generating measure instead 
of the integral of the latter because of the advantages of this choice to be discussed 
later.) Next, the LPM is obtained by recursively multiplying each rectangle with the 
generating measure k times. This is in complete analogy with the standard process of 
generating a multifractal, resulting in rn^^ rectangles, each associated with a linking 
probability Pij{k) equivalent to a product of k factors from the original generating pij 
given as 

k 

Pijik) = Ylpi.j,- (1) 

q=l 

In our convention k = 1 stands for the generating measure, thus, an LPM at = 1 is 
equivalent to the generating measure itself. The indices of the factors in ([T]) are given 
by 

z — 1 ) nr=i ° i^od m'^^'" 



k — q 



1, (2) 



where [a/h\ denotes the quotient (integer part) of a/6, the term nr=i omodm^"^ stands 
for subsequent calculation of the remainder after the division by m^"'', and an analogous 
formula can be written for the indices jq as well. (For g = 1, Eq.(l2]) simplifies to 
H — ~ l)/m^~^\ + 1). To obtain a network from the link probability measure, we 
distribute points independently, uniformly at random on the [0, 1] interval, and link 
each pair with a probability given by the Pij{k) at the given coordinates. (The above 
process is illustrated in Fig HI). 

The diversity of the linking probabilities Pij{k) (and correspondingly, the 
structuredness of the generated graph) is increasing with the number of iterations, just 
like in case of a standard multifractal. When considering the "thermodynamic limit" 
of this construction {k oo, — oo) we would like to keep the generated networks 
sparse, i.e., ensure that the average degree of the nodes, (d) remains constant. This can 
be achieved by an appropriate choice of the number of nodes as a function of k, using 
the following relation: 

{d) = N{k)Y,Y.P^A^)^^Ak). (3) 

i=i j=i 

where aij{k) denotes the area of the box i,j at iteration k. For simplicity, let us consider 
the special case of equal sized boxes aij{k) = nr"^^ . Due to the normalisation of the 
linking probabilities in this case the above expression simplifies to {d) = Nm~'^'', thus, 
to keep the average degree constant when increasing the number of iterations for a given 
generating measure, the number of nodes have to be increased exponentially with k. 

The above construction could be made more general by replacing the "standard" 
multifractal with the k-th tensorial product of a symmetric 2d function < W{x, y) < 




Figure 1. Schematic illustration of the multifractal graph generator, a) The 
construction of the link probability measure. We start from a symmetric generating 
measure on the unit square defined by a set of probabilities pij = pji associated with 
mxm rectangles (shown on the left). The generating measure is iterated by recursively 
multiplying each box with the generating measure itself as shown in the centre and on 
the right, b) Drawing linking probabilities from the obtained LPM. We assign random 
coordinates in the unit interval to the nodes in the graph and link each node pair /, J 
with a probability given by the probability measure at the corresponding coordinates, 
resulting in a graph drawn in the transparent green layer above the LPM. 

1 defined on the unit square. Although the resulting Wk{xi, ...,Xk,yi, ■■■,yk) = 
W{xi,yi) ■ ■ ■W{xk,yk) function is [0,1]^'' [0,1] instead of [0,1]^ -> [0,1], with the 
help of a measure preserving bijection between [0, 1] and [0, 1]*^ it could be used to 
generate random graphs in the same manner as with Pij{k). 

At this point we also note that omitting the normalisation condition ^^^Pij = 1 
for the generating measure gives the approach an additional flexibility which can come 
very handy in practical cases. Suppose that for a given setting of k, N and the LPM we 
would like to increase the average degree in the obtained graph beside preserving the 
relative ratios of the linking probabilities (expected degrees) of nodes falling into the 
different rows of the LPM. A very natural idea in this case is to multiply each element in 
Pij{k) with the same factor t] > 1, and use the resulting matrix for generating a random 
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graph in the same way. However, this multiphcative factor could be also introduced at 
the level of the generating measure instead, i.e., rj^^^Pij would also generate r]Pij{k) for 
the LPM. 



2.1. The degree distribution 

An important property of the MFNG is that nodes with coordinates falling into the same 
row (column) of the LPM are statistically identical. This means that e.g., the expected 
degree or clustering coefficient of the nodes in a given row is the same. Consequently, 
the distributions related to the topology are composed of sub-distributions associated 
with the individual rows. The degree distribution can be expressed as 

P^'\d) = Y,pf\d)k{kl (4) 

i=l 

where pf'\d) denotes the sub-distribution of the nodes in row z, and li{k) corresponds 
to the width of the row (giving the ratio of nodes in row i compared to the number of 
total nodes). These pf^ take the form of [21] 

pf)(rf) = i^^e"<'^'(^^)), (5) 

where {di{k)) denotes the average degree of nodes in row i. This is is given by 
{di{k)) = N{k)pi{k), where Pi{k) corresponds to the linking probability in row i, given 

by 

p.{k) = J2p^JikMk) (6) 

j 

In more general, if the linking probability at iteration k is described by Wk{x,y), then 
the expected degree of a node having a position x can be given as 

{dix))=N{k)wk{x), (7) 

where 

Wk{x) = / dyWk{x,y) (8) 



defines the Id projection of Wk{x,y) and is equivalent to the linking probability at 
position X. In case of the multifractal network generator this Wk{x) has a simple 
step-wise constant form, showing a step-wise surface getting rougher and rougher with 
increasing k. 



3. The isolated nodes and the multifractality of Wk{x) 

According to Sect J2.1[ . the degree distribution and the fraction of isolated nodes depend 
on the projection of the link probability measure Pij{k) to the x axis, (or equivalently, to 
the y axis), given by Wkix). It is known that almost all projections of a multifractal such 
as Pij{k) with an information dimension larger than 1 to a Id line result in measures with 
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Euclidean support [2S]- However, the projection to Wk{x) is unfortunately a special case, 
belonging to the minority of projections yielding the "almost all" instead of "every" in 
the previous statement. As we shall see shortly, Wk{x) is a multifractal itself with an 
information dimension smaller than 1, and this is the origin of the node isolation. 

In Figj2l we show Wk{x) for /c = 1 and /c = 2 in a setting with equal sized boxes 
in the generating measure. For simplicity we shall assume equal sized boxes in the rest 
of this Section. Since the linking probability inside each box is constant, the shape of 
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Figure 2. Projection of the LPM onto the edge of the unit square, resulting in a 
multifractal Wk{x) function. For simplicity wc assumed a 2 by 2 generating measure 
with equal box lengths. 

of Wk{x) is step- wise, consisting of intervals (corresponding to the columns of the 
LPM), and the height of step i is given by 



i=i 

However, the multiplicative nature of the construction is inherited by Wk{x) as well, 

k 



(9) 



q=l 



where 



"J' 



(10) 



stand for the heights of the steps at the generating measure (/c = 1), and iq is given by 
([2]). (This is demonstrated in Figl21 for k = 2). Thus, the evolution of Wk{x) with k is 
analogous to the standard construction of a multifractal embedded in the unit interval. 
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Note however that if 'Yl,ijPij = then Wk{x) is not normahsed, e.g., for equal box 
lengths j Wk{x)dx = m~'^^ . 

Multifractals are described by the q ordered generalised fractal dimension D{q) 
defined as follows (see e.g., Ref. [27]). Suppose that we divide the multifractal to boxes 
of size e, and the measure inside box i is given by pi. The function e) is defined as 

X(g,6)^5^p^ (12) 

i 

If e is varied, x(q') e) behaves as 

X(g,e)~e^('')(''-^). (13) 
Thus, D{q) can be given as 

D(q) = hm . (14) 

e->o [g - 1 Ine J 

In the special case of g = 1 we have zero in the denominator, thus we take the g — )■ 1 
limit and using I'Hospital's rule we obtain 

-^Pilnpi-Z^iMl/e). (15) 

i 

From the point of view of the degree distribution and the fraction of isolated nodes, the 
crucial q value is g = 1: When the number of iterations, /c — )■ oo, the fractal dimension 
of the support of the measure is given by D{q = 1). If it turns out that D{q = 1) < 1 
for the Wk{x) curve, then this means that the points giving the relevant contribution to 
the occurrence of links are concentrated on a fractal with a fractal dimension smaller 
than one, and thus, the majority of the nodes become isolated. 

For a multifractal on the unit interval defined by a self-similar multiplication 
process such as in case of Wk{x), (governed by Eqs.flQlfTTl)). the D{q) can be calculated 
analytically [27] as 



'^^^^= (g-l)L(l/m) ^" 



(16) 



=1 

For g = the above expression yields i5(g = 0) = 1, in agreement with the general 
picture of multifractals produced in a recursive multiplication process, where D{q = 0) 
equals the fractal dimension corresponding to a uniform generating object. For any 
multifractal in general, the D{q) values monotonically decrease with increasing g. Thus, 
at g = 1 in our case the D{q) can reach 1 only if D{q) = 1 for any g G [0, 1]. According 
to (fT6l) . this can be achieved only if r, = 1/m^ for all i. This means that unless the sum 
of probabilities in any row of the generating measure is the same, the D[q = 1) becomes 
smaller than 1, and the node isolation effect takes place. 



4. Rotated measures 



Our aim is to overcome the problem of the exact multifractality of Wk{x) by modifying 
the construction in such a way that the projection of Wk{x^y) determining the degree 
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distribution has Euclidean support. Meanwhile, we want to keep the LPM a highly 
variable function so that very different distributions (leading to very different kinds of 
networks) could be still achieved. A simple idea is to rotate the LPM with a given angle 
a as shown in Fig|3l, so that the direction of the projection determining the degree 
distribution no longer coincides with any of the special directions of the multifractal 
generation process. Thus, the construction of a random graph in this setting has the 
following main stages: we begin by generating a "standard" LPM for a chosen k, and 
"cut" a square rotated by an angle a from this measure as shown in Figl3l Since 
the diagonal of this newly introduced square does not coincide with the diagonal of the 
original LPM, we have to symmetrise the measure inside the rotated square with respect 
to its diagonal. Finally, we distribute points uniformly at random along both sides of 
the rotated square and link each pair with a probability found at the given coordinates 
in the "rotated coordinate system". (This last step is in complete analogy with the 
"standard" graph generation process). 




h 



Figure 3. A "standard" link-probability measure (oblique) with a rotated frame 
inside. The probabilities inside the rotated square have to be symmetrised along the 
diagonal. The linking probability of the green node can be calculated by summing up 
the probabilities "traversed" by the green line, each multiplied by the length of the 
intersection between the green line and the corresponding box. Due to the symmetry, 
this is equivalent to traversing first along the original line from the top until the 
diagonal is met, and then continue along the dashed line to the right. 



When examining the behaviour of this construction with increasing k, the number 
of nodes used in the graph generation is adjusted according to the {d) =const. criterion. 
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just as in case of the original settings. Due to the rotation and the symmetrisation, 
the polygons making up the link probability matrix are no longer arranged in a matrix 
like form, thus, ([3]) cannot be used straight away in this case. To calculate {d) we need 
to first introduce a unique indexing over the polygons inside the rotated frame, and 
evaluate 

{d) = Nj2Piai, (17) 

i 

where Pi and denote the probability and the area of the polygon i (triangles, 
quadrangles and pentagons) defining the rotated link probability measure. The size of 
the various distances and the equations for the most important lines needed to evaluate 
(fT7|) are given as a function of the rotation angle a in the Appendix. 

The expected degree of a node at distance h from the origin of the rotated square 
can be given by an expressions analogous to ([7]) as 

{d{h)) = N{k)wk{h), (18) 

where the linking probability Wk{h) has to be calculated along a line parallel to the side 
of the rotated square. By denoting the set of polygons intersected by this hne by ^kih), 
we can express Wk{h) as 

where Ij denotes the relative length of the intersections divided by the length of the 
side of the rotated frame (see Fig|3]). (These intersection lengths can be calculated with 
simple coordinate geometry, based on the equations given in the Appendix). Due to the 
symmetrisation of the rotated square, the boxes under the diagonal are rotated by 2a. 
Thus, in practise it is more simple to replace the summation ( fT9l) along a straight line 
by a summation along a broken line which is fully above the diagonal, as shown by the 
dashed green line in FiglHl 

The Wk{h) function is the analogue of the Wk{x) function for the rotated frame, 
and when a — )■ 0, Wk{h) — )■ Wk{x). As already noted in Sectj2l the Id projection of 
the original link probability measure, Wk{x), is always piecewise constant, where the 
constant intervals correspond to the columns of the link-probability measure. In case 
of Wk{h) the situation is a bit more complex. In FigHl we depict two close by node 
positions in the rotated square. The lines along which one has to calculate the linking 
probabilities intersect with the same boxes. Furthermore, the lengths of the intersections 
are the same for the two lines in most of the boxes, except for the sections marked by 
red. Thus, the linking probability of the two nodes given by Wk{hi) and Wk{h2) will 
be quite close to each other as well, with the difference coming from the few different 
intersection lengths. Now let us imagine that we fix one of these two node positions, 
and set the corresponding linking probability Wk{h) as a reference value. If we scan with 
the other node through a narrow interval of h values such that the corresponding line 
still intersects with the same boxes, then due to the linear change in the intersection 



Rotated multifractal network generator 



10 







_ PC 






















1 












/ 





























r 


























1 












i 



h] h2 



Figure 4. Two node positions in the rotated square, relatively close to each other 
shown by the green lines. They intersect with the same boxes, and in most cases the 
length of the intersections is the same as well. The different intersection lengths are 
marked by red. 



lengths, the change in Wk{h) with respect to the reference value will be linear as well as 
a function of the h difference between the two node positions. 

From this it follows that if we scan through the total range of possible h values, 
then the corresponding Wk{h) curve of the linking probabilities will be piece wise linear. 
The break points between the linear segments correspond to the h values where the line 
parallel to the rotated square boundary comes across a corner of the polygons defining 
the 2d link probability measure, and starts to intersect with a new polygon. Thus, 
to analytically calculate Wk{h) we need to evaluate f fT9l) only at the break points, and 
connect the results with linear segments. The h value of the break points (corresponding 
to polygon corners) can be most easily calculated by changing the coordinate system to 
be aligned with the rotated square, having an origin at the top left corner. The details of 
this coordinate transformation are given in the Appendix. In FigJS^ we check the above 
for a 2 by 2 generating measure at /c = 3 by calculating Wk{h) in both the break points 
and in two intermediate points between each adjacent break point pairs. Seemingly, the 
results for the intermediate points fall on the lines connecting the result for the break 
points. 

The degree distribution can be obtained from Wk{h) in three simple steps. The first 
step is the calculation of the distribution of the linking probability for the nodes, a^^\p). 
(By integrating a^^\p) as j^^ a^^\p)dp we receive the probability for a randomly chosen 
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node to obtain a linking probability between pi and p2, and the number of expected 
links on the node is given by the total number of nodes, A^, multiplied by its linking 
probability). The a^^\p) can be calculated from Wk{h) by a simple "projection" to 
the vertical axis as follows. Since we distribute the nodes uniformly at random along 
the side of the rotated square, the ratio of nodes falling into an interval of [/ii,/i2] is 
simply — hi)/b, where b denotes the length of the side of the rotated square. Thus, 
a linear segment of Wk{h), stretching from hi to /i2 contributes to a^^\p) with a "step" 
ranging from pi = Wkihi) to p2 = Wk{h2) with a height of (/12 — hi)/ {b\p2 — pi\). (For 
an illustration see FigjSb). 

The second step in the calculation of the degree distribution is the transformation of 
(T^^\p) into the distribution of the expected degrees for the nodes, p^^\d). The difference 
between the degree distribution and p^^\d) can be illustrated in the graph generation 
process: the expected degree of a node is simply its linking probability given by (fT9l) 
multiplied by iV, however, since the links are drawn randomly, its actual degree may 
become smaller or larger than that at the end of the link generation process. The p^^^ (d) 
can be obtained from a^^\p) by a simple "stretching" in the horizontal direction, i.e., 
for any pi and p2 

a'^^\p)dp= / p'^^\x)dx, (20) 

pi Jxi=Npi 

where the integral on the right hand side corresponds to the probability for a randomly 
chosen node to have an expected degree falling between di = Npi and ^2 = ^^2- We 
note that in case of the original settings without any rotation, both a^'^\p) and p^'^^d) 
are given by trains of delta spikes with varying weights. 



{p)= Y,m5,,r>m^ (21) 
1=1 

i=l 

where Pi{k) denotes the linking probability in column i of the original LPM given by 
([6]). In contrast, for the rotated measures both a^'^^p) and p^'^^d) take a step-wise form 
instead of delta spikes, (see e.g., FiglSb). 

The final step is to transform p^^\d) into the degree distribution. Since the links 
are drawn independently of each other in both the original and the rotated settings, this 
can be achieved by taking the convolution of p^^\d) with a Poisson-distribution as 

pW(rf) = j dxp^''\x)'^e-\ (23) 

Based on the method detailed above, in the next Section we compare the evolution of 
the degree distribution with k in the original settings and in the rotated scenario, (with 
a special focus on the ratio of the isolated nodes). 
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Figure 5. a) Checking the Wk{h) of a 2 by 2 generating measure at fc = 3. The 
continuous hne shows the piece wise hnear hnking probability as a function of h 
obtained by connecting the values obtained at adjacent break points with straight 
lines. For each section we also calculated the linking probability in two intermediate 
points as well, the results are shown by the symbols, b) By projecting Wk{h) to the 
vertical axis one obtains the distribution of the linking probabilities, crij)), which is 
step- wise. (Note that in order to emphasise the connection to Wkih), we have plotted 
p on the vertical axis and cr(p) on the horizontal axis), c) The distribution of the 
expected degrees, p'-'^''{d) is obtained by inflating a{p) according to the (d) = Np 
relation. We have also plotted the corresponding degree distribution, p'^^\d)^ with 
dashed lines. 



5. Applications 

5.1. Degree distribution 

For the comparison between the original settings and the rotated scenario we chose 
the 2 by 2 generating measure shown in Figj6|L. with a starting Nk=i = 30, since the 
ratio of isolated nodes becomes significant quite fast without the rotation of the LPM 
in this case. The generating measure and A'^^^i fixed the average degree of the nodes 
to (d) = 7.5. From (d) we can calculate the number of nodes at the higher k values in 
the original setting from ([3]). In case of the rotated scenario we use the same (d), and 
calculate the number of nodes for any k value from ( |T7|) . 

In FigJBJ) we show the distribution of the expected degrees for the nodes, p^''\d), on 
logarithmic scale for the original settings, obtained from (l22|l . whereas FiglHb-d show the 
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Figure 6. a) A 2 by 2 generating measure with equal box lengths, b) The distribution 
of the expected degrees, p'^''^ (d), (black, solid lines), and the degree distribution, p^''^ (d), 
(red, dashed lines) in the original settings for fc = 1, 2, 4, 6, 8 on logarithmic scale. Since 
p^''^ (d) corresponds to a delta-spike train according to (P^ . we used logarithmic binning 
for dealing with the singularities, c) The p'^^\d) (black, solid lines), and p'^^^d), (red, 
dashed lines) distributions at k = 1,2,4,6,8 for a rotated frame at rotation angle 
a = 0.5 degrees, d) The p^^^d) (black, sohd lines), and p^^\d), (red, dashed lines) 
distributions at fc = 1, 2, 4, 6, 8 for a rotated frame at rotation angle a ~ 22.5 degrees. 



same function for rotated frames at rotation angles a = 0.5 degrees and a = 22.5 degrees, 
respectively. We also plotted the corresponding degree distributions with dashed lines, 
however, the difference between the original and the rotated scenario is much more 
salient for p'^^\d). As explained in SectJH, p^'^^d) consists of delta-spikes in case of 
the original settings (we have used binning in order to plot this singular function), and 
according to FigJHb, as k is increased, the distribution gets wider, and a significant 
part of it is shifted under d = 1. This means that as we increase k, for larger and 
larger part of the nodes the expected degree becomes smaller than one, thus, the node 
isolation effect takes place. In contrast, Figsj6t-d show a different behaviour. Although 
the distributions become wider with increasing k here as well, this tendency is much 
less pronounced compared to the a = case. Furthermore, in case of Fig|6li the major 
part of the distribution stays above d = 1 for the examined k values. 

In FigJT^ we show the ratio of isolated nodes as a function of the number of 
iterations. Due to the reduced spreading in the degree distribution with k the rapid 
increasing tendency of p{d = 0) present in the original settings is modified to a very 
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slowly increasing tendency for the rotated scenarios. The p(d = 0) at iteration = 10 is 
displayed in FiglTb as a function of the rotation angle a, showing an overall "U" shape 
with a minimum around 22.5 degrees. As pointed out previously, when a — ?■ 0, we 
recover the original settings of the MFNG. Although there seem to be many pairs of a 
values with very similar p{d = 0) values, we note that each a defines a different setting 
with unique attributes (e.g., the area of the rotated square affecting the average degree 
and the family of the Wk{h) curves are always unique for each a). The overall properties 
of a rotated setting (degree distribution, ratio of isolated nodes, etc.) change smoothly 
with the rotation angle for any a > 0. However, as argued in SectJH theoretically and 
shall be examined in Sect J5.2[ numerically, the behaviour of the D{q) curve does show 
a drastic change when switching from the a = original setting to an a > finite 
rotation angle. This change from a multifractal D{q) to a non- multifractal one is not 
expected to be affected by any commensurability effect in rotation angles. 




1 2 3 4 5 6 7 8 9 10 5 10 15 20 25 30 35 40 45 

k « 



Figure 7. a) The ratio of isolated nodes as a function of the number of iterations for 
the examined generating measure (shown in FigE^). b) The ratio of isolated nodes at 
fc = 10 as a function of the rotation angle a. 



One of the big advantages of the original MFNG was that it provides a flexible tool 
for generating random graphs with realistic properties. Although examining to what 
extent this feature is affected by the rotation of the LPM is beyond the goals of this 
article, we present an example where power-law like degree distribution was achieved in 
the rotated scenario in the Appendix. 

5.2. Measuring the D{q) curve 

The results shown in Sect J5.ll are very promising, since the rotation of the LPM 
drastically reduced the ratio of the isolated nodes in the studied example, especially 
in the case of larger rotation angles. However, an even more reassuring way for checking 
the effect of the modification of the MFNG is the measurement of the D{q) curves 
corresponding to the Wk{h) functions. As described in Sect|3l, from the D{q) at g = 1 
we can deduce the behaviour of the fraction of isolated nodes, i.e., if D(g = 1) = 1 in 
the /c — )■ oo limit, than the isolated nodes cannot become dominant. 
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Before proceeding to the D{q) curves, we have one important note. In practise 
we are always deahng with multifractals at a finite number of iterations, which have a 
lower size bound (lower length bound in our case) depending on the number of iterations. 
Below this lower size bound they are not any more structured and, thus, this size bound 
provides a lower bound for the size e of the boxes with which we cover the multifractal 
when measuring D{q) as described in SectJS] 

In case of the Wk{x) of the original setting (without any rotation) this lower bound 
in e is given simply by the width of the rows (columns). In case of the Wk{h) of the 
rotated measure the situation is a bit more complex. First of all, when compared to the 
Wk{x) of the original setting at the same number of iterations, Wk{h) can contain much 
smaller segments. However, there are many adjacent segments with the same or almost 
the same slope of the linking probability, which can be united into single large segment, 
as shown in FigEl Thus, before the application of the D{q) measuring procedure above, 

2e-19 
1.8e-19 
1.6e-19 
1.4e-19 
1.2e-19 

le-19 

8e-20 

6e-20 

0.1 0.101 0.102 0.103 0.104 0.105 

h 

Figure 8. A small part of the Wk{h) of the rotated link probability measure at fc = 8, 
a = 10.5 degrees. The points correspond to the boundaries of the segments in which 
Wkih) is strictly linear. However, for many adjacent segments the slope is actually the 
same or nearly the same. 

these segment joining were carried out, and the lower bound of e was set to the average 
of the length of the resulting segments. This yielded still a much smaller value than in 
case of the original measure. 

Interestingly, for the rotated frames when plotting lnx(g,e) as a function of Ine, 
the curves seem to consist of two subsequent linear segments with different slopes. In 
contrast, the lnx(g,e) obtained from Wk{x) in the rotation free settings shows a linear 
behaviour as a function of Ine. The two types of xiQ^ ^) curves are shown in FigjQl The 
D[q) is obtained from the slope of these curves, which is straight forward in case of the 
Wk{x) of the original settings. However, which part of the lnx(g,e) curve should we fit 
in case of the rotated measures? According to the definition given in ( !T4|) D{q) should 
be evaluated in the e — )■ limit, thus, we use the slope of the lower part of the \nx{q, e) 
curves in the rotated scenario for measuring D{q). 

In FigJTOl we show the results obtained for the D{q) curves at various rotation 
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Figure 9. a) The x{li^) obtained for the Wk{x) of the original setting at fc = 8 and 
g = 1.5 on logarithmic scale (black symbols), together with a power- law fit (red line), 
b) The x(g,e) of the Wk{h) of the rotated frame at rotation angle a = 10.5 degrees 
with the same parameters fc = 8 and (j = 1.5 on logarithmic scale (black symbols). 
In contrast to the original settings, when fitting ^(g, e) with a power-law, the lower 
part and the upper part of the function yield different exponents, as shown by the 
continuous red line and the dashed green line respectively. 



angles. (For a comparison, the D{q) curve of the original settings is shown as well). 
From these figures it seems that D{q) takes a trivial "non- multifractal" form at already 




Figure 10. a) The D{q) curve of the Wk{x) function, corresponding to the original, 
"rotation-free" settings, obtained from ([T5)) . The D{q) curves measured in the rotated 
scenario for fc = 4, 6, 8 at rotation angles a = 0.5°, a = 10.5° and a = 22.5° are shown 
in panels b), c), and d), respectively. 
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k = 8, and it is very close to 1 at g = 1. The "precision" of the numerical D(q) 
determination process at g = 1 was tested on other Id functions where we know that 
D[q = 1) = 1. According to that, for the rotated measures the numerically obtained 
D{q = 1) values are equal to D{q = 1) = 1 within error bound for k = 8. 

6. Summary 

In summary, we investigated the node isolation effect of the MFNG from a new point 
of view. It had become clear that this phenomena is very closely related to the 
multifractality of the 1 dimensional projection of the LPM determining the degree 
distribution. According to general theorems concerning multifractals, the projection 
in question is a particular one, and in contrast, the vast majority of the other Id 
projections of the LPM do not bear multifractal properties. Based on this observation 
we introduced a slight variation of the original MFNG method, involving the rotation of 
the LPM with a given angle. Due to the rotation, the projection determining the degree 
distribution is no longer a special projection in any aspects, thus, the node isolation 
effect is expected to disappear due to the lack of multifractality. The empirical studies 
support the theoretical reasoning above: the order q generalised fractal dimension D[q) 
in case of the projection of the LPM related to the degree distribution became trivial, 
and for the numerically accessible range of the number of iterations the fraction of 
isolated nodes showed a drastic reducement when compared to the original settings 
without any rotation. 
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Appendix 

Al. Lengths 

For the length of the interval a in Fig {11] we can write 



a 

= tana, 

1 — a 



(24) 



thus, 



tana 



(25) 



1 + tan a 



Similarly, we can write for b 



a 

- = sma, 





(26) 
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t(x) 




Figure 11. The geometry and the definition of the various lengths and sizes in the 
rotated hnk probabiUty measure. 



thus, 



Furthermore, 



whereas 



tan a 



1 



sina(l+tana) cos a + sin a 



c = a tan a 



tan^ a 
1 + tan a ' 



sin(a + 7r/4) sma + cosa 

a = atanfa + 7r/4) = a -— = a . 

cos(a + 7r/4) cosa — sina 



(27) 

(28) 
(29) 



A2. Equations of the different lines 

In the coordinate system shown in FiglTTl the equation of the top side of the rotated 
frame can be written as 

t(x) = 1 + c — xtan(a). (30) 

Similarly, the equation of a parallel line below at a distance of h is given by 

r(x) = 1 + c — /icos(a) — /isin(a) tan(a) — xtan(a). (31) 
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The equation of the diagonal can be expressed as 

g{x) = 1 + d - xtan{a + tt/A). (32) 
Finally, the equation of the left side of the rotated square is given by 

s{x) = a + xtan{n/2 — a), (33) 

and the equation of a parallel line shifted by h to the right (corresponding to the dashed 
green line in FigJTTl) can be written as 

r (x) = a + [x — h cos a — h sin a tan a) tan(7r /2 — a). (34) 

The equations given above can be used to calculate where a given line intersects with 
a given box boundary of the original link-probability measure, and from the location of 
the intersection points we can calculate the lengths/areas of the intersections between 
the lines and the boxes. 

To calculate (d) we have to calculate the area of each intersected box in the top 
triangular part of the rotated square. The boundary lines of this triangle can intersect 
with the boxes as summarised in Fi^T2l: 

• in case the box intersects with only one line, then this divides it into either two 
trapezoids (FigJT2k). or into a triangle and a pentagon (FigJT2b). 

• in case the box intersects with two lines, then there are still only three intersection 
points (instead of four), since the boundary lines must intersect each other on the 
boundary of the box. The remaining two unshared intersection points can take 
qualitatively three different positions: they can both fall on a side adjacent to 
the side of the shared intersection point (FigJT2b). they can both fall on the side 
opposite to the shared intersection point (FigJT2li). they can fall on two opposite 
sides adjacent to the side of the shared intersection (FigJT2b). or they can fall on 
two adjacent sides (FigJT2f). 



A3. Changing coordinate system 

According to FigJT3l. the new coordinates of a general point x, y can be given as 
X — a 

X = --— + (t(x) — y) sin(a), (35) 

cos(a) 

y = (t(x) - y) cos(a), (36) 

where t{x) denotes the equation of the topside of the rotated square in the standard 
coordinate system, given by ( l30i) . 

AJi^. Generating skewed p^^\d) in the rotated scenario 

In Figim we show an example where power-law like degree distribution is generated 
in the rotated scenario. According to the plots, a slight rotation preserves the degree 
distribution almost completely, whereas for larger rotations the skewed nature of p*^*^-' {d) 
is slowly disappearing. 




Figure 12. The six type of intersections of the boxes with the boundary hnes of the 
top triangular part of the rotated square. 
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